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The hole-doped antiferromagnetic spin-i two-leg ladder is an important model system for the 
high-Tc superconductors based on cuprates. Using the technique of self-similar continuous unitary 
transformations we derive effective Hamiltonians for the charge motion in these ladders. The key 
advantage of this technique is that it provides effective models explicitly in the thermodynamic limit. 
A real space restriction of the generator of the transformation allows us to explore the experimentally 
relevant parameter space. From the effective Hamiltonians we calculate the dispersions for single 
holes. Further calculations will enable the calculation of the interaction of two holes so that a handle 
of Cooper pair formation is within reach. 
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I. INTRODUCTION 

The hole-doped antiferromagnetic spin-i two- leg lad- 
der is used as a model system, theoretically as well as ex- 
perimentally, for the two-dimensional cuprate supercon- 
ductors. Because the spin ladder is quasi one-dimensional 
it can more easily be treated by numerical and analytical 
approaches than the full two-dimensional model. For this 
reason the spin ladder was subject to many theoretical 
investigations in recent years^— . From a strong cou- 
pling perspective it is reasonable to describe the doped 
spin ladder by a t- J-mode l 17 ' 18 . 

In this paper, our aim is to derive an effective Hamil- 
tonian which describes not only the motion of the mag- 
netic degrees of freedom, i.e., the triplonsi£, but also of 
the charges, i.e., of the doped holes. This task is more 
challenging than the description of the triplons alone be- 
cause the two kinds of excitations interact strongly and 
their energy ranges are not separated but they strongly 
overlap. Hence the excitations are not infinitely long- 
lived but they decay, at least in certain regions of the 
Brillouin zone. 

The technique employed here is the approach of self- 
similar continuous unitary transformations (SCUT). In 
this approach the mapping of the true ground state to a 
vacuum of elementary excitations is constructed system- 
atically, details will be explained in Sect. IIIII The basic 
idea is to implement a transformation that adjusts itself 
during the procedure of the diagonalisation depending on 
the current form of the Hamiltonian at this instant of the 
continuous transformation. The change of the Hamilto- 
nian induced by the transformation is determined by the 
current magnitude of the non-diagonal elements. The 
flowing Hamiltonian is given in second quantization, i.e., 
it has a certain structure in terms of elementary creation 
and annihilation operators 2 ^. This formulation incorpo- 
rates the linked cluster property 2 ^ automatically. The 
SCUT systematically defines a set of differential equa- 
tions for the prefactors of the monomials of the creation 
and annihilation operators. 



n + 1 



J\\ , til 



Jj_, t± 

— 1 



FIG. 1: Spin-i ladder consisting of the upper leg u, the lower 
leg 1 and the rungs n. The magnetic couplings are the nearest 
neighbor coupling on the rungs J± , the nearest neighbor cou- 
pling along the legs Jii and the four spin coupling Jb for the 
ring exchange. The nearest neighbor hopping is described by 
the constants t± (on the rungs) and tu (along the legs). The 
symmetry axis of the ladder is depicted by the dashed line. 



The key issue is to construct the flow to the effective 
Hamiltonian in a robust way. Due to the energetically 
overlapping states this goal required modifications^ 2 - with 
respect to the previously performed unitary transforma- 
tions. Thereby we achieved the successful calculation of 
the dispersion of a single hole. 

In Sects. IIVI and [V] we investigate the dispersions of 
single hole excitations obtained from the effective Hamil- 
tonians, which are results of the SCUT. We also compare 
our results to the results of other methods. 

Sects . IVTl and I VIll conclude this article with a summary 
and perspectives for future work. 



II. SPIN-i LADDERS 



The 



ladder is realised in the so called tele- 



spin- | 

phone number compounds (Sr, La, Ca, Y) 14 Cu2404i— . 
For these systems superconductivity can be detected un- 
der high pressure 2 ^. Interactions between the ladders are 
weak because they result from 90° exchang o 25 i 26 . Thus 
the ladders can be considered to be isolated from one 
other. Even if the interladder coupling is taken into ac- 
count, the strong frustration of the lattice causes the sys- 
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tern to be effectively one-dimensional2£. The dispersions 
of the magnetic excitations of the complete layer with 
interladder coupling are similar to those of a single spin 
ladder. Important magnetic properties of the ladders are 
investigated in Refs. [ll- lill0lll4 

The focus of the present article is the lightly doped 
spin ladder because this will allow us to address 
the dynamics of charges. Of course, this system 
has been investigated previously by many techniques, 
for instance exact diagonalizatior&i^, density-matrix 
renormalizatio n 6 i 7 i 15 , high-order series expansion^, di- 
agrammatic approache a 9 i 12 ' 13 , and quantum Monte 
Carloii. Hence much is known about the dynamics of 
a single hole, i.e., its dispersion, its spectral density and 
to some extent also about the bound states of two holes. 
The aim of the present paper is to establish an approach 
which generates an effective model for the motion of holes 
and their interaction with magnetic excitations and other 
holes. Such a model can be expressed in terms of creation 
and annihilation operators which create and annihilate 
holes, i.e., charge excitations, and triplons, i.e., magnetic 
excitations. The techniques used so far do not yield such 
an effective model but numerical data for eigen energies 
and spectral weights. For this reason we here use the 
complementary approach of self-similar continuous uni- 
tary transformations. 

The spin sites in our model (see Fig.[TJ are the 3d x 2 -y 2 
orbitals of the copper atoms coupled via the 2p x or 
2p y orbitals of the oxygen atoms, which hybridise with 
the copper 3d x 2 _ y 2 orbitals so that an antiferromagnetic 
superexchange^ is possible. The Hamiltonian in the 
usual form of a i-J-model, for derivations see Ref. 
and references therein, reads 



H = 



^ Sn tU ■ S n> l + J|| >5n,Q ' 5Vi+l,a (1) 

n n,a— u.l 

n 



^^1- ^ ] c a-,a,n c a,a,n ~ ^|| ^ ] C <x,a,n C er,a,n±l ' 
n,a— u,l n,a— u.l 

The nearest neighbor terms with the couplings J± and 
J|| (see Fig. [TJ are not sufficient to describe the magnetic 
interactions in the system. Because the hybridization 
path around the CU4O4 square plaquettes is strong, the 
influence of the four-spin interactions belonging to these 
plaquettes (i.e. two neighboring rungs) is not negligible 30 . 
These interactions are referred to as ring exchange (also 
cyclic exchange) with the coupling Jjj (see Fig. [T]) . Ac- 
tually the complete ring exchange also includes two-spin 
terms for all two-spin combinations of the four spins on 
two neighboring rungs. Yet the contributions from the 
terms coupling the spins along the rungs and parallel 
to the legs are merged with the nearest neighbor terms 



so that the coupling constants J± and Jy include these 
contributions, whereas the terms coupling the spins di- 
agonally can be neglected because their prefactor is only 
of the order of 3% of J±— . The complete representation 
of the ring exchange using spin operators can be found 
in Ref. [32- 

The hopping of the spins is described by the constants 
ij_ and i|| (see Fig. [l}. The creation and annihilation op- 
erators and c already incorporate hardcore properties 
so that double occupancy is forbidden. 

We define the dimensionless quantities 

X = X D = — , Aj_ = — , A|| = (2) 

J± J±_ Jj_ Jj_ 

which are normalized to J±. 

It is convenient to treat the ladder as a one-dimensional 
chain of rungs. There are nine possible local states on a 
rung n of the ladder. 
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The states without holes are the singlet (|3aj) and the 
triplet states (|3bll3d|) . for which we use the representa- 
tion from Ref. |3J because the symmetry in spin space is 
manifest in this representation. The states consisting of 
one hole and one spin i (|3ell3hp have the quantum num- 
bers r, which is the parity of the state with respect to 
the symmetry axis of the ladder (see dashed line in Fig. 
[1]), and a, which indicates whether the S z component of 
the spin is up or down. The double hole state ([31]) is 
neglected here because only one-hole states will be con- 
sidered. The local energy of \2h) is also larger than the 
energy of the other states. Therefore \2h) is expected to 
be also less important for ladders with low, but macro- 
scopic doping, though it will influence the quantitative 
results for hole-hole scattering and hole-hole binding. 

For the application of the continuous unitary transfor- 
mation a representation using the creation and annihi- 
lation operators for the local states (|3all3h[) is indicated. 
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The complete Hamiltonian (QJ in the representation us- 
ing these operators can be found in Ref.[3J. Note that the 
local rung states incorporate hardcore properties. Also 
the fermionic states exhibit additional hardcore proper- 
ties which simply result from the fact that each rung can 
only be in one particular state. 

The excitations are actually dominated by the lo- 
cal states as long as the couplings between the rungs 
are small. The magnetic excitations are triplons, i.e., 
triplet states dressed with the magnetic interactions with 
their environment^ 9 -. We do not want to use the term 
"magnon" for these magnetic excitations because this 
term is usually associated with quasiparticles in systems 
that exhibit long-range magnetic order (which is not the 
case for the spin ladder). Moreover triplons feature the 
appropriate threefold degeneracy based on their 5=1 
character—. 



III. SCUT 

We use the method of continuous unitary transforma- 
tions (CUT) to obtain an effective Hamiltonian 3 ^— . In- 
stead of applying only a single constant unitary trans- 
formation that diagonalises the Hamiltonian at once or 
several constant unitary transformations successively, a 
unitary transformation depending on a continuous pa- 
rameter is applied to the Hamiltonian. This transfor- 
mation adjusts itself permanently during its application. 
Discrete transformations must be known explicitly before 
we can apply them, whereas for the continuous transfor- 
mation it is sufficient to set up the infinitesimal generator 
of the transformation. The unitary operator U depending 
on the continuous parameter I transforms the Hamilto- 
nian by 

H(l) = U{l)HU\l). (4) 
The generator of the transformation is defined by 

vd) - ^fV(0- (5) 

The generator defines the properties of the transforma- 
tion. The derivative of H with respect to / is given by 
the so-called flow equation 



dH(l) 
dl 



[V(1),H(1)} 



(6) 



This is actually a system of differential equations (gener- 
ically highly coupled) for the coefficients of the operators 
appearing in H{1). To achieve the effective Hamiltonian 
the flow equation has to be solved. For a self-similar 
application of the CU T 22 ' 38 ' 39 the Hamiltonian is repre- 
sented by a sum of different operators di , which are prod- 
ucts of local operators in second quantization, multiplied 
by prefactors gi. 



The operators di determine the structure of H(l). This 
is the reason why this formulation of the CUTs is called 
self-similar. In particular, we express the operators 
di in terms of creation and annihilation operators of 
charge excitations (holes) and spin excitations (triplons) 
in the manner introduced previously 2 ^. A reference 
state, the vacuum, has to be chosen with respect to 
which these creation and annihilation processes are de- 
fined. In our case, this state is the product of local 
singlets on the rungs of the ladder. For the undoped 
ladder this state is a convenient starting point as seen 
in many high-order perturbative^2r— and diagrammatic 
approaches 12,13 . Since we focus in this work on single 
hole dynamics the use of the product of local singlets on 
the rungs as reference vacuum is well justified. 

The differential equations for the are given by the 
flow equation in the form 



%(0 

dl 



3,k 



(8) 



(7) 



Because an infinite system yields an infinite number of 
differential equations which cannot be solved numerically, 
we have to apply a truncation to H(l) so that a closed 
system of differential equations is obtained. We choose 
a real space truncation because the correlation length is 
finite for the spin ladder due to the energy gapi^. The 
concrete truncation scheme is explained in Sect. IIII CI 

It should also be noted that the truncation is only 
restricted to operators. There is no truncation of the 
Hilbert space of the states. Thus although operators af- 
fecting higher numbers of particles are omitted, the num- 
ber of particles that can be treated is arbitrarily large. To 
illustrate this we consider the action of a term in second 
quantisation, e.g., a\a^a kl which annihilates one particle 
and creates two. This term acts not only on the one- 
particle subspace, but also on the subspaces with more 
than one particle, e.g. it changes a four-particle into a 
five-particle state. 

The solution of the flow equation is a rather straight- 
forward numerical integration. The convergence is mon- 
itored during the integration. Since all contributions to 
the generator decrease in case of convergence, the abso- 
lute values of the concerning coefficients are squared and 
summed. The square root of this sum is defined as resid- 
ual off-diagonality (ROD). It is the norm of the generator 
and also a measure for the convergence. The ROD is ex- 
pected to tend to zero for I — > oo. Note that the term 
"off-diagonality" is meant in a broader sense, i.e. the def- 
inition of the generator determines which elements shall 
be kept for I — ¥ oo and these elements are defined as "di- 
agonal" parts of the Hamiltonian. The RODs depicted 
in the following are always normalized to the initial ROD 
at I = 0. 

With the decrease of the ROD the wanted effective 
Hamiltonian is approached. If the relative ROD falls be- 
low a certain threshold specifying the precision of the 
result (usually 10~ 8 ), the integration can be consid- 
ered to be completed. The number of coefficients treated 
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I = finite I I = oo 



FIG. 2: (Color online) Schematic example for a transforma- 
tion of H(l) induced by the pc generator: Each coloured (n c , 
n a )-block represents contributions to H creating n c particles 
after annihilating n a particles. 



in our actual calculations is of the order up to 10 5 for the 
doped spin ladder. 

A. Generator 

The particle conserving or pc generator 

V P c,i,j(l) = signfe - qj)H itj (l) (9) 

was introduced by Mielke for band matrices 4 ^ for q t = i 
and simultaneously by Knetter and Uhrig for many-body 
Hamiltonians where qi counts the number of elementary 
excitation o 44 ' 45 . The indices i and j label the transition 
from state j to state i. These states are eigenstates of 
an operator Q, which counts the number of elementary 
excitations, i.e., quasi-particles. The eigen value of Q for 
the state i is 

The differences between the pc and the Wegner gen- 
erator are discussed by means of simple Hamiltonians in 
Ref. liH In particular the convergence behavior is exam- 
ined. The Wegner generator ?7wcgncr = [Ho,Hl^ always 
leads to a fixed point, but it does not necessarily aim at 
a quasi-particle picture because degeneracies may hinder 
the diagonalisation. The pc generator is not sensitive to 
degeneracies, but the induced transformation does not 
always converge for infinite systems. 

Due to the sign function the generator @ contains only 
terms of the Hamiltonian that change the particle num- 
ber. Terms which decrease the particle number acquire 
an additional minus sign. The resulting transformation 
yields an effective Hamiltonian for I — > oo which con- 
serves the particle number (see Fig. [2]). The asymptotic 
behavio r 22 ' 43 of the non-particle-conserving terms, i.e., 
<7i qji is dominated by 

« Mum,/, q^UL.-h (10) 

Because H iy j(l) tends to zero for I — > oo if qi ^ qj, the 
transformation tries to sort the eigenenergies ffy accord- 
ing to the particle number i.e., 

signfe - qj )(Hu(l) - H jtj (l)) > (11) 

holds true for large I. 
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FIG. 3: (Color online) Schematic example for an overlap be- 
tween one- and two-particle energies. The black line is the 
one-particle dispersion, the solid cyan (gray) lines are the 
boundaries of the two-particle continuum and the dashed lines 
are the energies of the two-particle states which are actually 
affected by the terms an a l + An a ™ arm a m a n+AnCin- The cyan 
(gray) dashed lines do not pose a problem because they do not 
cross the one-particle dispersion. Convergence problems can 
only be induced by the terms that act on the states with the 
cyan (gray) dotted dispersions, which cross the one-particle 
dispersion. 



Overlapping energies of states with different quasi- 
particles can hinder the convergence of the flow. This 
problem can be solved by excluding terms from the 
generator—. The ground state or gs generator 

%M,j(0 = ( S 3,o - &i,o)Hi,j(l) (12) 

introduced in Ref. [22| includes only terms that couple to 
the vacuum \i = 0) of excitations so that in the effective 
Hamiltonian only this vacuum is decoupled from the re- 
maining Hilbert space. Then it has become the ground 
state. This is sufficient for the calculation of the hole 
dispersion for the spin ladder if the ground state can be 
continuously linked to the state without magnetic exci- 
tations for a given hole configuration. This is the case for 
a single hole at a given momentum close to the minimum 
of its band. 



B. Real Space Restriction of the Generator 

But even the gs generator encounters convergence 
problems if the interaction between the rungs becomes 
too strong, i.e., x = 1, \± = Ay > 2. Hence we have de- 
veloped a generator adaption with a real space restriction 
of the generator terms. 

Many terms that remain after the real space truncation 
of the Hamiltonian act on states with discrete energies 
within the continua (see Fig. [3]). The divergence of the 
SCUT in case of an overlap is not induced by all terms 
that mediate between overlapping continua, but only by a 
part of these terms. Thus a stricter real space truncation 
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of the Hamiltonian can induce convergence. However, a 
stricter truncation causes a larger error for the results of 
the SCUT. 

Therefore we do not apply a stricter truncation to the 
Hamiltonian but we choose a stricter generator. We em- 
phasize that a restriction of the generator does not imply 
an approximation. It only changes the direction of the 
unitary rotation. This restriction of the generator is not 
based on the total extension of the terms, but on the ex- 
tension of the creation operators and on the extension of 
the annihilation operators separately. 

Let us consider the example of an overlap between one- 
and two-particle energies shown in Fig. |3l The terms 

a n a n+An a m an d a m a n+An a n are responsible for transi- 
tions between the one- and the two-particle subspace. We 
choose to restrict the generator based on |An|, which is 
the distance between the two particles which are created 
or annihilated. Because the one-particle state does not 
have an extension, it is not relevant for the restriction 
of the generator. All terms that have a larger |An| than 
a certain An max are excluded from the generator. Note 
that they are still part of the Hamiltonian as long as they 
meet the truncation criteria for the Hamiltonian. 

The transformation induced by this restricted genera- 
tor does not try to sort all eigenenergies (cf. Eq. [11]) , but 
only those which are captured by the terms in the gen- 
erator. Therefore the flow may also converge in case of 
overlapping energies. The price to be paid is that the sub- 
spaces affected by the omitted terms are not completely 
decoupled from the remaining Hilbert space. Thus cither 
an additional diagonalization or a diagrammatic analy- 
sis has to be applied to the effective Hamiltonian from 
this SCUT or the results have to be considered as upper 
limit for the actual results. The restriction of the gener- 
ator can be applied to the pc generator or to any of its 
adaptions, e.g., to the gs generator. 



C. Truncation Scheme 

The finite energy gap A of the triplonsi^ is equivalent 
to a correlation function which is exponentially decreas- 
ing with respect to the distance. Hence a truncation in 
real space is appropriate. The extension in real space 
shall be used as a measure for the physical importance of 
a term of the Hamiltonian. In our quasi one-dimensional 
spin ladder the extension of a term is defined as the dif- 
ference between the smallest and the largest rung index 
of the local operators within the term. 

The simplest way of truncating would be one maximal 
extension in real space for all terms. Terms exceeding 
this limit would be omitted. But this approach is not 
reasonable in our case. The number of possible terms in- 
creases much more strongly with the maximal extension 
for terms consisting of more local operators. However, 
terms consisting of less operators are usually more im- 
portant. Hence higher extensions should be taken into 
account for terms with less operators. For instance, the 



coefficient of a one-particle hopping term consisting of 
two local operators is usually larger than the coefficient 
of a two-particle interaction term consisting of four local 
operators if they both have the same extension. There- 
fore, for the pure triplon terms different maximal exten- 
sions d n are defined in units of the rung distance where n 
is the total number of local operators of the term under 
study. Because only terms changing the particle number 
by two or conserving the particle number occur in the 
initial magnetic Hamiltonian, only the d n with even n 
are meaningful^. 

By h n i we denote the maximal extension for n' hole 
state operators. It does not matter if the term concerned 
contains additional triplon operators. For terms with 
triplon and hole state operators, h n > affects only the hole 
state operators. The total maximal extension t n n for 
these mixed terms depends on the number of triplon op- 
erators n" . The parameter N t denotes the maximal num- 
ber of triplons interacting with holes. Note that t n » with 
odd n" have to be taken into account because the triplon 
number is changed by an odd number if the hole state 
parity is altered. 

Additionally a maximal triplon number N and a max- 
imal hole number Nh are defined for the operators, i.e., 
terms affecting higher triplon numbers than N or higher 
hole numbers than Nh are completely omitted. This 
truncation scheme is not compulsory. Other schemes 
could be implemented which use different classifications 
for the groups of terms sharing the same maximal exten- 
sion. 

Since the size of the system of differential equations 
for the doped ladder grows exponentially with increasing 
extensions, the truncation for our calculations is fairly 
strict. The parameters N = 4, c?2 = 10, c?4 = 6, d^ = 4 
and g?8 — 3 were used for the pure triplon terms because 
these maximal extensions are sufficient for the undoped 
case up to x — 1. 

For the terms including hole operators the truncation 
used for the SCUT is given by N h = 1, N t = 2, h 2 = 3, 
t\ = 6, ta = 6, <3 = 5, ti = 5, t§ = 4 and te = 4. 
All results shown in the sequel are obtained with these 
parameters unless stated otherwise. 



IV. HOLE DISPERSIONS WITHOUT RING 
EXCHANGE 

A. Calculations with Unrestricted Generators 

The dispersion of a single hole in the absence of triplons 
can be easily derived from the effective Hamiltonian by 
means of Fourier transformation. The part of the Hamil- 
tonian to be diagonalised 

h 2 

#lh = Y. hd a l<?,n+d a r,a,n (13) 

d= — h 2 T,a,n 
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FIG. 4: (Color online) One-hole dispersion for x — xa — 
calculated with the pc generator; left: A = 0.25, right: A = 
0.5. The curves coincide well with the series expansion results 
from 
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FIG. 5: (Color online) One-hole dispersion for x = xu = 0, 
A = 1 calculated with the gs generator and the pc generator. 
The results for the even bands coincide so that the curves 
cannot be distinguished. The cyan (gray) shaded area is the 
continuum formed by one triplon and one even hole state. 



contains the one-hole terms restricted by /12 and charac- 
terised by the coefficients b^. The one-hole dispersion 

hi 

Wlh,r(fc) = ^0 + X! 2&dCOS(dfc) (14) 
d=l 

only depends on the parity r and it is degenerate con- 
cerning the spin a. 

If A11 and Aj^ are small while x = x\j = 0, the dispersion 
reads 

" lh -g (fc) = ^TA ± + A||Cos(fc)+0(Ai,Af,A ± A || ). (15) 

In the following we only consider the isotropic case A^ = 
An = A. For small values of the parameters x, and A 
the deviations from Eq. (|15p are actually small. In Fig.0] 
two values of A are considered for x = x = 0. The result 
for A = 0.25 already exhibits deviations from the relation 
(|15|l for small hopping constants. Both bandwidths are 
smaller than 2A = 0.5. The odd band is narrower than 
the even band. Moreover the odd band is shifted upwards 
by less than A, while the even band is shifted downwards 
by slightly more than A. 

The increase of A to 0.5 yields a result with obvious 
deviations relative to Eq. (IT51) . The odd band becomes 
lower and narrower with growing A in this region. The 
cosine shape of both dispersions is - for A = 0.25 as well 
as for A = 0.5 - not deformed by higher harmonics. The 
results of the series expansion- exhibit the same behavior 
for these parameters in good agreement with the SCUT 
results. However, it has to be pointed out that the con- 
vergence of the SCUT is much worse for A = 0.5 than 
for A = 0.25. The residual off-diagonality (ROD) defined 
in Sect. IIIII and used as a measure for the convergence 
is decreasing very slowly for A = 0.5. While the ROD is 
smaller than 10~ 6 at U± = 200 for A = 0.25, it is still 
larger than 10~ 3 at U ± = 200 for A = 0.5. Both RODs 



are decreasing exponentially for large I. Theses RODs 
are not shown here. 

For x = xjj — 0, A = 1 the pc generator yields results 
for the even band that still agree very well with the series 
expansion. For the odd band there are no series expan- 
sion results available. Fig. [3] shows the comparison of 
the results for the pc generator and for the gs generator. 
The even band results lie on top of each other, but the 
odd band exhibits deviations for k > 0.577T. While the 
gs generator produces an almost featureless dispersion 
(compared to the even band), the pc generator causes a 
pronounced maximum at k — ~k. These differences are 
due to the position of the lower boundary of the contin- 
uum formed by one triplon and one even hole stated (see 
also Fig. [3]). This is the odd continuum due to the odd 
parity of the triplon. An overlap between the odd dis- 
persion and this continuum is present for the pc result. 
The deviating odd band from the gs generator avoids this 
overlap. For both generators the even and the odd band 
cross at k sa 0.467T. While the even band keeps its cosine 
shape, the second harmonic for the odd band is no longer 
negligible for both the pc and the gs generator. 

For x = x\j — 0, A = 1 the transformation does not 
converge for the pc generator, while the gs generator in- 
duces convergence (see Fig. [5]). The overlap hinders the 
convergence for the pc generator. The kink at I J± ps 120 
in the gs ROD is probably due to numerical inaccuracies 
that are amplified via a feedback within the flow equa- 
tion. Since the relative ROD is already w 10~ 7 at the 
kink we can stop the transformation at this point and 
neglect the remaining off-diagonal terms, i.e., we con- 
sider the transformation to be converged. For the calcu- 
lation of the hole dispersions the pc SCUT was stopped 
at IJ± ps 50 where the ROD is still ~ 10 -2 and thus not 
negligible. Therefore the result from the gs generator 
is more trustable than the result from the pc generator. 
The parameters x = x\j = 0, A = 2 lead to divergence 
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FIG. 6: ROD for the SCUT for x = xa = 0, A = 1 induced by 
the pc generator (solid line) and by the gs generator (dashed 
line). 



for either generator. 



1. Finite Coupling Ratio x = J\\/J± 



FIG. 7: (Color online) Left panel: One-hole dispersion for 
x — 0.5, xa = 0, A = 0.25 calculated with the pc generator. 
The series expansion result for the odd band is also shown, 
while the result for the even band is not available. Right 
panel: One-hole dispersion for x = 0.5, xn = 0, A = 0.5 
calculated with the pc generator and the gs generator. For 
case a) in the legend the usual truncation was used. For case 
b) in the legend hi was increased to 4, t\ and t% to 8 as well 
as tz and fci to 6. The dotted curves are extrapolations from 
the series expansion data. 



Here we want to consider x > 0. For A = the energy 
of the hole states, which is independent of momentum 
k and parity r, decreases with increasing x. Already 
for small finite A the deviations from the simple cosine 
shape appear for the odd band. This can be seen in Fig. 
[TJwhere the one-hole dispersions for x = 0.5 and A = 0.25 
are depicted. The odd band agrees well with the series 
expansion; only slight deviations at k — occur. For the 
even band no series expansion data are available. 

For x = 0.5 and A = 0.5 , see also Fig. the result for 
the even band agrees again well with the series expansion 
result^, but the odd band behaves differently. Only for 
large k the behavior is similar although also in this region 
the SCUT result is slightly lower. The series expansion 
result exhibits a local maximum at k = and a global 
minimum at k 0.587T, while the SCUT result hardly 
changes in the region < k < J. The gs generator 
yields the same result as the pc generator apart from 
minimal deviations (< 2%) at k = 0. However, the gs 
generator allows us to extend the truncation scheme: hi 
was increased to 4, t\ and ti to 8 as well as £3 and ti to 6. 
In the result the shape of the odd band changes mainly 
for small k. Thus the odd dispersion still changes with 
increasing maximal extensions. 

We can understand these differences between SCUT 
and series expansion by looking at the position of the 
continuum formed by one triplon and one even hole state 
(see Fig. [HJ . The series expansion result crosses the lower 
boundary of this continuum at k rs 0.34-7T. Assuming that 
the series expansion yields a result close to the actual 
odd hole dispersion it is cogent that the SCUT based 
on the pc generator is not able to sort the eigen values 
properly for small k because creating one triplon does 
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FIG. 8: (Color online) Comparison of the odd hole dispersion 
and the continuum formed by one triplon and one even hole 
state for x = 0.5, xn = 0, A = 0.5. The odd hole dispersion 
was derived from SCUT using the pc generator (cyan/gray 
solid) and by series expansion (cyan/gray dotted). The cyan 
(gray) shaded area is the continuum (derived from the same 
SCUT calculation). 



not necessarily increase the energy. 

Since the problematic overlap concerns the zero- and 
the one-triplon space, the use of the gs generator is no 
remedy here. This may come as a surprise at first since 
it was shown in Ref. [22] that the gs generator always im- 
plies a robust flow. This statement applies if the hole 
is put in its band minimum because this is by construc- 
tion the lowest eigen state of the ladder doped with one 
hole. But for other energies of the hole dispersion this 
does not need to be the case so that the energy of one 
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FIG. 9: (Color online) One-hole dispersion for x = 1, xa = 0, 
A = 0.5 calculated with the pc generator and the gs gener- 
ator. The dashed curves are extrapolations from the series 
expansion data. 



FIG. 11: (Color online) One-hole dispersion for x = 1, xn = 0, 
A = 1 calculated with the pc generator and the gs generator. 
The dashed curve is the result from the series expansion, for 
which only the even band is available. The dashed-dotted 
curves are exact diagonalization results. 
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FIG. 10: (Color online) Comparison of the odd hole dispersion 
and the continuum formed by one triplon and one even hole 
state for a; = 1, xn = 0, A = 0.5. The odd hole dispersion was 
derived from SCUT using the pc generator (cyan/gray solid), 
using the gs generator (cyan/gray dotted) and by series ex- 
pansion (cyan/gray dashed). The cyan (gray) shaded area is 
the continuum (derived from the same pc SCUT calculation). 



hole may lie higher than the energy of one hole and one 
triplon. Hence the gs generator yields mainly the same 
result and suffers from the same problems as the pc gen- 
erator. A robust extrapolation of the series can still yield 
reasonable results because the series is not fully sensitive 
to the overlap of states at higher values of the expansion 
parameters. 

For x = 1 and A = 0.5 (see Fig. [9]) the SCUT results 
for the odd band again show distinct deviations from the 
series expansion results. The even band exhibits a max- 
imal deviation at k = of only « 1%. The deviations of 
the energetically higher odd band are larger, but still well 
tolerable for the gs generator. All in all the deviations 
are not as pronounced as for x = 0.5 and A = 0.25. A 
comparison with the continuum formed by one even hole 



and one triplon (see Fig. I10[) shows that an overlap ex- 
ists around k = it for the series expansion result and for 
the gs result. This explains again the deviations in this 
region. As this overlap is not as strong as the overlap for 
x = 0.5 and A = 0.5 the deviations are smaller. Even if 
there is no actual overlap, the continuum is at least very 
close. Hence it is to be expected that the actual odd hole 
dispersion is lowered for k w n due to this reason. 

For x = 1 and A = 1 (see Fig. QT]) no one-hole disper- 
sion with odd parity is available from series expansion. 
Thus we look at the even band including results from 
exact diagonalization. For the exact diagonalization a fi- 
nite ladder with 14 rungs was examined and the resulting 
eigenenergies were fitted by a sum of three cosine terms 
for both bands, see Eq. (fT4|) . The pc result for the even 
band differs completely from the gs result. For the odd 
band the pc result lies below the gs result and the devia- 
tions grow with increasing momentum. The comparison 
with the series expansion and with the exact diagonaliza- 
tion suggests that the gs result is more reliable because 
the deviations are smaller than 1% for the even band. 
Also for the odd band the gs result is closer to the exact 
diagonalization result. The odd band is more difficult to 
assess due to the vicinity to the odd continuum. 

It should be noted that the diagrammatic approach 
from Ref. HH yields bands that show qualitative devia- 
tions from the SCUT and from the exact diagonalization 
concerning the shape for small k. The authors state that 
this regime actually exceeds the applicability of their ap- 
proach. The quantum Monte Carlo result from Ref. [ll] 
is in very good agreement with our result. For the even 
band the agreement is even excellent and comparable to 
the agreement between SCUT and exact diagonalization. 

The odd band result from the gs SCUT enters the con- 
tinuum formed by one triplon and one even hole state for 
k > 0.837T (see Fig. H"2"j) . while the pc result lies always 
below this continuum. This overlap is also present in the 
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FIG. 12: (Color online) Comparison of the hole dispersions 
and the continua formed by one triplon and one hole state for 
x = 1, asp = 0, A = 1. Left: The odd hole dispersion was 
derived from SCUT using the pc generator (cyan/gray solid) 
and using the gs generator (cyan/gray dotted). The cyan 
(gray) shaded area depicts the continuum with one even hole 
and one triplon derived from the same gs SCUT calculation. 
Right: The even hole dispersion was derived from SCUT using 
the pc generator (black solid), using the gs generator (black 
dotted) and by series expansion (black dashed). The black 
shaded area depicts the continuum with one odd hole and 
one triplon derived from the same gs SCUT calculation. 



exact diagonalization result and in the quantum Monte 
Carlo result. A comparison of the even one-hole dis- 
persion and the even continuum formed by one triplon 
and one odd hole state (see also Fig. [12]) supports the 
assumption that the gs result is more reliable than the 
pc result. The dispersion induced by the gs SCUT ex- 
hibits a shape that appears as if it were formed by the 
lower boundary of the approaching continuum. This is 
plausible because the vicinity of the continuum lowers the 
dispersion. The pc result, however, stays away from the 
continuum at k — and at the boundary of the Brillouin 
zone. The overlap with the continuum around k w 0.447T 
occurs in the pc calculation only due to truncation. 

At this point we compare the convergence behavior for 
x = 0.5, A = 0.5; x = 1, A = 0.5 and x = 1, A = 1 for 
the pc generator. The corresponding RODs are depicted 
in Fig. 1131 All three curves exhibit a kink with non- 
converging behavior afterwards. This is typical for cu- 
mulating rounding errors because of a symmetry breaking 
due to numerical inaccuracies. Such a symmetry break- 
ing is not unlikely because the full spin symmetry could 
not be used explicitly. 

It is surprising that the ROD for x = 1, A = 0.5 reaches 
a lower minimum than the one for the lower parameter 
set x — 0.5, A = 0.5. This seems to contradict the ex- 
pectation that larger values of the interdimer processes 
imply a more difficult CUT. But inspecting the slope of 
the RODs in Fig. Q2] before reaching their minima one 
sees that these slopes, i.e., the convergence velocities, are 
indeed largest for the smallest values for the interdimer 



FIG. 13: ROD for the SCUT calculated with the pc generator 
for various parameters. 




0.4 0.6 
k[jt] 

FIG. 14: (Color online) One-hole dispersion for x = 1, ajp = 0, 
A = 2 calculated with the pc generator and the gs generator. 
The graph also shows the exact diagonalization data (ED) 
and the series expansion data (SE)£. For the latter only the 
even band is available. The even band results are black, the 
odd band results are cyan (gray). 



processes. 



2. Strong Hopping A = t/Jx > 1 



For x = 1 and A > 1 the parameters are entering a 
region which is expected to reflect realistic relations of 
the constants in the telephone number compounds. The 
results for x = 1 and A = 2 are shown in Fig. Q31 Again 
we compare with results from series expansion and from 
exact diagonalization for 14 rungs. 

Approximate analytic results were obtained in Ref. 
by perturbation theory improved by a variational ansatz. 
These results lie even above the series expansion results 
but confirm the qualitative shape for the even band. The 
pc result for both the even and the odd band is again 
very distinct from the other results comparable to x = 1 , 
A = 1. The gs result, however, is in better agreement 
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FIG. 15: (Color online) Comparison of the hole dispersions 
and the continua formed by one triplon and one hole state 
for x = 1, xu = 0, A = 2. Left: The odd hole disper- 
sion was derived from pc SCUT (cyan/gray solid), gs SCUT 
(cyan/gray dashed) and by exact diagonalization (cyan/gray 
dotted). The cyan (gray) shaded area depicts the continuum 
with one even hole and one triplon (derived from the same 
gs SCUT calculation). The cyan (gray) dash dot dot line is 
the lower boundary of the continuum formed by two triplons 
and one odd hole. Right: The even hole dispersion was de- 
rived from pc SCUT (black solid), gs SCUT (black dashed), 
by exact diagonalization (black dotted) and by series expan- 
sion (black dashed-dotted) . The black shaded area depicts 
the continuum with one odd hole and one triplon (derived 
from the same gs SCUT calculation). The black dash dot dot 
line is the lower boundary of the continuum formed by two 
triplons and one even hole. 



with the data from the series expansion and especially 
with the exact diagonalization result in accordance with 
our previous observations. 

Apart from the pc result the dispersions exhibit the 
same features. The even band has a global maximum 
at k = and a local maximum at k = n, while it is 
vice versa for the odd band. Because both bands lie in 
the same energy range, they cross in the middle between 
k = and k = tt. The exact diagonalization predicts 
the crossing to be at k s» 0.487T, but the gs SCUT finds 
the crossing at k s» 0.557T. The even band calculated by 
series expansion is located above both the series expan- 
sion and the gs SCUT result for all k. This is further 
evidence that the extrapolation used to correct the bare 
series underestimates the lowering of the band induced 
by the hybridisation with the hole-triplon continuum. 

Let us consider the continua formed by one hole and 
one triplon. The continua consisting of one hole and one 
triplon are compared to the one-hole dispersions in Fig. 
1151 The continua do not overlap with the hole disper- 
sions, but they are very close to them. The only excep- 
tion is the series expansion result for the odd band which 
actually exhibits an overlap with the continuum formed 
by a triplon and an even hole state. 

The behavior of the ROD yields further evidence why 
the gs results should be preferred to the pc results for this 



FIG. 16: ROD for the SCUT for x = 1, x D = 0, A = 2 in- 
duced by the pc generator (solid line) and by the gs generator 
(dashed line). 



parameter regime. Fig. [16] depicts the evolution of the 
ROD during the flow for both generators. Not only that 
the pc SCUT converges extremely slowly, the shape of the 
curve for IJ± < 100 is an indicator for a problem with 
respect to the sorting of the eigenenergies. It is a generic 
behavior of the ROD that the sorting of the eigen values 
(cf. Eq. (jTTJ)) is reflected by features for small values of I. 
When the sorting is completed the ROD decreases expo- 
nentially with a constant rate. For the gs generator the 
decrease of the ROD attains this rate not later than at 
I J± = 5. Before this point the decrease is slower—. The 
kink of the gs ROD at IJ± « 110 with the subsequent 
upturn is again most probably due to cumulating numer- 
ical inaccuracies. This is no real problem because the gs 
ROD has already fallen below a value of less than 10~ 16 
at IJ± ss 110 and can hence be neglected. However, the 
pc ROD exhibits several humps and a pronounced rise 
at IJ± — 34 before a decrease with a constant rate is 
achieved. This is a typical indication for a suppressed di- 
vergence that would actually occur for a less strict trun- 
cation. If such a behavior occurs the transformation is 
susceptible to truncation errors. The physical origin of 
these problems is the strong overlap between the one- 
hole-one-triplon continuua and the one-hole-two-triplon 
continua (see Fig. IT5"|). 

B. Calculations with Restricted Generators 

Results from exact diagonalization are available for 
x = 1 and A = 3, but even the gs generator does not 
induce convergence for this case. However, if we apply 
the generator restriction introduced in Sect. IIIIBI to the 
gs generator, it yields convergence for An max ^ 2. Then 
some coupling between one hole and one hole plus one 
triplon lh — » lh + It is not eliminated completely. Be- 
cause the Hamiltonian is not diagonalised with respect 
to the terms omitted from the generator, the hole disper- 
sions we obtain from a Fourier transformation are only 
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FIG. 17: (Color online) Comparison of the one-hole disper- 
sions from the restricted generator with the gs and the exact 
diagonalization results. Left panel: x = 1, xn = 0, A = 2. 
Right panel: x = 1, xa — 0, A = 3. The black lines represent 
the even hole dispersion and the cyan (gray) lines represent 
the odd hole dispersion. 



upper limits for the actual result. But this is a minor 
problem since the remaining matrix elements are small. 
Another way of improvement (not followed here) would 
be to use self-consistent Born approximation to deal with 
the remaining coupling lh — > lh + It. 

Here we compare the results for upper limits from the 
restricted generator for x = 1 and A = 2 with the gs 
results before we investigate the results for x = 1 and 
A = 3. The left panel of Fig. [T7] shows the comparison 
of the one- hole dispersions for x = 1 and A = 2. We 
see that for An max *S 2 the upper boundary from the 
restricted gs generator is close to the result from the full 
gs generator. In the right panel of Fig. [T7] the results from 
the restricted gs generator are compared to the exact 
diagonalization results for x = 1 and A = 3. For the even 
hole dispersion the agreement between the result from 
the restricted gs generator and the exact diagonalization 
result is almost perfect. Also the agreement for the odd 
hole dispersion is goodi^i The deviations are comparable 
to the deviations of the result by the full gs generator 
from the exact diagonalization result for x = 1 and A = 2. 

The investigation of the ROD shows that the restricted 
gs generator yields a faster convergence than the full gs 
generator for x = 1 and A = 2 (see Fig. [TS]) . For x — 
1 and A = 3 the ROD diverges for the full generator, 
while the restricted generator induces convergence (see 
also Fig. [IH]). Note that the kinks of the RODs in Fig. 
[IS] with an increase afterwards are again probably due to 
a feedback of numerical inaccuracies. But all the kinks 
appear at values where the ROD is already smaller than 
10~ 12 . So the flow can be considered to be converged at 
the kinks for practical purposes. 

The terms left out from the restricted generator still 
contribute to the Hamiltonian after the transformation. 
These contributions yield an estimate of the difference 
between the actual energy and the upper boundary 
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FIG. 18: (Color online) ROD for the SCUT calculated with 
the full gs generator and with the restricted gs generator for 
different parameters. 




FIG. 19: (Color online) Comparison of gs SCUT (solid lines) 
and exact diagonalization for various finite ladders with L 
rungs (discrete points) for x = 1, A = 1. The lower boundaries 
of the continua calculated by gs SCUT are shown as dashed 
lines. The results for the even hole state are black, the results 
for the odd hole state are cyan (gray). 



for the energy resulting from the restricted generator. 



For the 432 terms of the form at 

a T,a,m^a,n a ; ., 

cients is « 0.48 J^; its square root yields 



or 



n+An * ne sum over their squared coeffi 



0.16J X . 

To understand the deviations between SCUT and exact 
diagonalization we also compare our results to the com- 
plete spectrum from exact diagonalization. In Fig. [19] 
the results of the exact diagonalization for various finite 
ladders with L rungs are compared to the results from 
the gs generator for x = 1 and A = 1. For the even band 
the lowest lying eigen value can be clearly distinguished 
from the larger eigen values, which are the precursor of 
the continuum. This holds true for all momenta. In the 
case of the odd band the lowest eigen value is very close 
to the higher ones for large momenta. Our calculation of 
the continuum predicts that the dispersion merges with 
the continuum in this region. Also the quantum Monte 
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FIG. 20: (Color online) Comparison of gs SCUT (solid lines) 
and exact diagonalization for various finite ladders with L 
rungs (discrete points) for x — 1, A = 2. The lower boundaries 
of the continua calculated by gs SCUT are shown as dashed 
lines. The results for the even hole state are black, the results 
for the odd hole state are cyan (gray). 




FIG. 21: (Color online) Comparison of SCUT from the re- 
stricted gs generator with An max = 2 (solid lines) and exact 
diagonalization for various numbers or rungs L (symbols) for 
x — 1, A = 3. The dashed lines are estimates for the lower 
boundaries of the continua based on the result from the re- 
stricted gs generator. The results for the even hole state are 
black, the results for the odd hole state are cyan (gray). 



Carlo result for the spectral weight^ exhibits no peak 
below the continuum around k ~ 7T. 

Our result for the even band is in excellent agreement 
with the exact diagonalization result. Around k rs 0.257T 
our result for the odd band lies above the exact diagonal- 
ization result while it lies below the exact diagonalization 
result around k rs 0.67T. In Fig. [20] we compare the re- 
sults of the exact diagonalization with the results from 
the gs generator for x = 1 and A = 2. For the even band 
for large momenta the lowest lying eigen value is clearly 
distinguishable from the larger eigen values, which are 
again the precursor of the continuum. The same is true 
for the odd band for small momenta. But for small mo- 
menta in case of the even band and for large momenta in 
case of the odd band the lowest eigen value is very close 
to the higher ones. Hence the distinction between hole 
dispersion and continuum becomes difficult in these re- 
gions. It is even questionable if they are actually distinct 
or if the dispersion merges with the continuum. In these 
regions the deviations between SCUT and exact diago- 
nalization are the largest. We conclude that the SCUT 
suffers from truncation errors if the dispersion runs close 
to continua or even enters them. Also the lower boundary 
of the continuum that we calculated from the gs result is 
higher at k — for the even band and at k — tt for the 
odd band than we would expect from the exact diagonal- 
ization data. 

For x = 1 and A = 3 the exact diagonalization results 
are shown in Fig. [21] for various finite ladders with L 
rungs. The odd hole dispersion is more difficult to dis- 
tinguish from the continuum than for x = 1 and A = 2. 
Also for small momenta the lowest eigen value is very 
close to the higher ones. The hole dispersions from the 
restricted gs generator with An max = 2 is also depicted 
in Fig. [2TJ In the region 0.57T < k < 0.87T this dispersion 
lies below the exact diagonalization data. 
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FIG. 22: (Color online) One-hole dispersions with even parity 
(black) and odd parity (cyan/gray) for x — 0, X± — Ay = 0.5 
and various values for xq. 



V. HOLE DISPERSIONS WITH RING 
EXCHANGE 

Since the ring exchange is needed for an adequate 
description of the experimentally available systems, we 
also investigate the influence of the ring exchange on the 
doped ladder, see Ref. [l4| and references therein. A typ- 
ical value is x\j rj 0. 2 14 i 47 . As for the examination of 
the anisotropic hopping we take the reliable results for 
x = 0, A ± = A|| = 0.5; x = 0.5, A_l = A|| = 0.25 and 
x — 1, Aj^ — A || = 2 without ring exchange as starting 
point for our investigation. 

For the first case x = 0, A^ = A|| = 0.5 the influence 
of the ring exchange is the weakest (see Fig. l22j) . The 
odd band hardly changes. It is slightly lowered - only 
the maximum is decreased more strongly. For the even 
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FIG. 23: (Color online) One-hole dispersions with even parity 
(black) and odd parity (cyan/gray) for x = 0.5, Ax = Ay = 
0.25 and various values for xu- 
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FIG. 24: (Color online) One-hole dispersions with even parity 
(black) and odd parity (cyan/gray) for x = 1, A_l = Am = 2 
and various values for xrj. 



band, however, the changes are more pronounced. It 
is strongly lowered so that a band crossing occurs and 
the shape changes. A local miminum appears at k = 
developing into the global minimum with growing xrj. 

Fig. [23] depicts the case x = 0.5, A^ = Am = 0.25. 
Both bands are lowered stronger for small k than for 
large k. Also the shape of both bands changes. While 
the minimum of the even band moves to k « 0.657T for 
&□ = 0.2 and a local maximum occurs at k — tt, the local 
maximum at k = of the odd band moves to fc ~ 0.557T 
and a local minimum at occurs k = 0. 

For X± = A|| =2 and x — 1 only the gs results are 
discussed because the pc generator yields no conclusive 
results for these parameters without ring exchange. The 
resulting one-hole dispersions are shown in Fig. (3D The 
increase of xrj yields a lowered dispersion for both bands. 
The shape of the bands is conserved, but for the odd 
band the decrease is pronounced around k = tt. Even an 
increase of the energy can be observed for the even band 
around k = tt. 



VI. SUMMARY 

The aim of the present paper was to establish an 
approach for doped low-dimensional quantum antifcrro- 
magnets which yields an effective model for the motion 
of charges. To this end, we employed a self-similar con- 
tinuous unitary transformation (SCUT) to lightly doped 
spin ladders. The method is particularly suited to pro- 
vide effective models beyond numerical results for special 
quantities. Spin ladders are a particularly well suited 
system because their magnetic excitations are very well 
understoo d 14 ^ 6 . 

We are able to calculate the one-hole dispersions by 
means of SCUT. The agreement with results from series 
expansion is very good for small parameters. Even in 
the regime x = 1, X\j = 0, 2 < A < 3 the agreement 
with the exact diagonalisation results is still very good 
for the even band and satisfactory for the odd one. The 
hole dispersions are strongly influenced by the triplons. 
Because the one-hole dispersion with odd parity has a 
larger local energy which is influenced more strongly by 
the higher lying continuum, it changes more explicitly 
with increasing hopping constants. Therefore the devia- 
tions from the simple cosine shape are more pronounced 
for the odd band. These deviations grow if either the 
magnetic coupling or the strength of the hopping is in- 
creased: 

• The broadening of the odd band is slowing down, 
then turning into a narrowing. Finally the shape 
changes completely under the influence of the sec- 
ond harmonic so that the maximum at k = is 
only local, the total maximum occurs at k = tt and 
the minimum lies between and tt. 

• The shift upwards is also slowing down and then 
turning into a shift downwards due to the influence 
of the continuum formed by one triplon and one 
even hole state. In the regime around x = 1 and 
X± = A|| = 2 the average odd dispersion is approx- 
imately as low as the average even dispersion. 

For the even band the deviations from the cosine shaped 
dispersion consist essentially in the growth of the sec- 
ond harmonic. For x = 1 and A^ = A|| =2 the mini- 
mum moves from k = tt into the region k w A local 
maximum at k = tt occurs. But the absolute maximum 
remains at k = 0. 

The combination of these effects for x = 1 and Aj^ = 
Ay = 2 yields a crossing of the two dispersions. The cross- 
ing point lies between the minima of the bands. How- 
ever, the pc generator is not applicable for the SCUT 
if the parameters are in this region. This is suggested 
by the deviations from the series expansion and from ex- 
act diagonalisation as well as by the peculiar convergence 
behavior of the SCUT. 

In the regime x ~ 1 and Aj_ = A|| > 1 the pc generator 
is no longer suitable because the convergence of the flow is 
hindered by the overlap between the one-hole-one-triplon 
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continuua and the one-hole-two-triplon continua. The 
pc results deviate from the series expansion and exact 
diagonalisation results. Furthermore the convergence is 
very slow and exhibits features that indicate problems in 
the sorting of the eigen values. 

The remedy is to use the gs generator which only 
decouples the zero-triplon subspace from the remaining 
Hilbert space. Then the hole dispersions are similar to 
the exact diagonalisation results. Although there are 
small deviations from the series expansion, the agreement 
is astonishingly good. Because the exact diagonalisation 
studies a ladder with fourteen rungs, finite size effects are 
present so that a part of the deviations stem from them. 
The remaining deviations are probably due to truncation 
errors . Another aspect in favour of the gs generator is 
the satisfactory convergence behavior. 

The treatment of the case Aj_ = Aii > 2 is problem- 
atic even for the gs generator because the convergence is 
hindered by the overlap between the odd one-hole disper- 
sion and the continuum formed by one even hole and one 
triplon. To achieve convergence in this regime we have 
developed the following restriction for the gs generator. 
A term affecting the hole-triplon continuum is omitted, 
if the distance between the triplon and the hole state on 
which the term acts is larger than An max . Strictly, the 
effective Hamiltonian yields an upper boundary for the 
hole dispersions. The comparison of the results from the 
full gs generator and from the restricted gs generator for 
x = 1 and A^ = Am = 2 shows that the upper bound- 
ary given by the result from the restricted generator is 
already very close to the result from the full generator 
for An max = 2. For x — 1 and Aj_ = A|| =3 this re- 
striction induces convergence, while the flow diverges for 
An max > 2. The estimates we obtain for the hole dis- 
persions agree again well with the exact diagonalisation 
results. 

If x — l,X]j = and A = A^ = Aii = 3 or 
A = A^ = A|| =2 the bands exhibit a similar shape 
and relative position to each other so that still a crossing 
at k Ri 0.57T occurs. But the energy is lowered and the 
bandwith of both bands is increased by a factor of ps 1.5 if 
A is changed from 2 to 3. The results from the restricted 
gs generator yield an estimate for the even band which is 
in good agreement with the exact diagonalisation result, 
while the odd hole dispersion exhibits deviations that are 
very similar to the deviations between the exact diago- 
nalisation and the full gs generator for A^ = A|| = 2. 

The case x = X± = \\\ = 0.5 is a special one. Both the 
pc and the gs generator (even with increased maximal 
extensions) exhibit deviations from the series expansion 
results for the odd hole state. These deviations stem 
apparently from the closeness (or even overlap) of the 
continuum formed by one even hole and one triplon. If 
the lowering of the odd band is overestimated by the 
SCUT or underestimated by the series expansion is not 
clear. 

The ring exchange leads to a lowering of the hole dis- 
persions. This effect is least pronounced for the even 



band around k = ir. The deformation of the band shape 
is most pronounced for the odd band. 

Summarizing, the essential achievement of this paper is 
the development of a unitary transformation which allows 
us to disentangle the motion of the doped charges, the 
holes, and the magnetic excitations, triplons. Thereby 
an effective model for the motion of the holes and the 
triplons has been derived. The challenge was to system- 
atically derive this effective model even in the experi- 
mentally relevant regime A ~ 3 where the hopping takes 
about three times the value of the magnetic exchange 
constants. This challenge could only be met by modifi- 
cation of the unitary transformation in the spirit of what 
was done in Ref . [22J. 



VII. OUTLOOK 

Ladders doped with more than one hole have not yet 
been treated by the SCUT. To do so the approach pre- 
sented here has to be extended to include also interaction 
terms between two holes. This step will be subject of fu- 
ture research because it will address the key question how 
large the attractive forces between two doped charges are. 
In view of the Cooper pair formation in doped cuprate su- 
perconductors this is a very interesting issue. Of course, 
numerical and diagrammatic results for bound states of 
two holes exist&^ i 13 ' 15 . But so far no effective model has 
been derived which incorporates the attractive potential 
explicitly. 

The fact that we could realize a systematic mapping 
of the initial doped spin ladder to an effective model for 
the motion of single excitations is encouraging for the 
next step incorporating the hole-hole interaction. The 
additionally required resources, for instance the increase 
in the number of coefficients, are significant. We estimate 
that about six times more terms have to be kept track 
of. Still this should be realizable in the near future. 

We expect that a substantial gain in the understanding 
of Cooper pair formation and thus of superconductivity 
in doped Mott insulators will be possible. 
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FIG. 25: (Color online) Finite size scaling of the exact diag- 
onalization results for x = 1 A = 1 in case of the odd band 
at k = and at k — 0.5n. The cyan (gray) circles depict 
the exact diagonalization results for finite ladders with var- 
ious numbers of rungs L. The solid black line is a linear 
extrapolation with respect to j- and the dashed line is an ex- 
trapolation based on exponential saturation (see Eq. (jAlfl ~) . 
For comparison the result from the gs SCUT is also shown 
(solid cyan/gray). 




FIG. 26: (Color online) Finite size scaling of the exact diag- 
onalization results for x = 1 A = 2 in case of the odd band 
at k = and at k — 0.5n. The cyan (gray) circles depict 
the exact diagonalization results for finite ladders with var- 
ious numbers of rungs L. The solid black line is a linear 
extrapolation with respect to j- and the dashed line is an ex- 
trapolation based on exponential saturation (see Eq. (lAHl ). 
For comparison the result from the gs SCUT is also shown 
(solid cyan/gray). 



Appendix A: Finite Size Scaling of the Exact 
Diagonalisation Results 

To understand the deviations between SCUT and exact 
diagonalization in case of the odd band we also study the 
finite size scaling of the exact diagonalization. First we 
consider the exact diagonalization results for x = 1 and 
A = 1 shown in Fig. [19] Around k ps 0.257T our result for 
the odd band lies above the exact diagonalization result, 
while it lies below the exact diagonalization result around 
k w 0.6tt. 

An investigation of the finite size scaling for the exact 
diagonalization data in these regions is difficult because 
we have two points at most for an extrapolation. Thus 
we consider k = 0, where an extrapolation is conclusive. 
The result of this extrapolation is used to support an 
extrapolation in the region where the deviations are ob- 
served. The left panel of Fig. shows the finite size 
scaling for k — using two kinds of extrapolation. The 
first is a simple linear extrapolation with respect to 
while the second assumes an exponential saturation with 
increasing L so that 

Awae"T (Al) 

holds true for the difference Acj from the limit for L — » 
co. The correlation length £ is determined to be « 4.32 
by this extrapolation. The correlation length can be used 
to apply the second extrapolation also for k = 0.5tt where 
only two results are obtained by exact diagonalization. 

The finite size scaling for k = 0.5tt in case of the odd 
band is investigated by both extrapolations in the right 
panel of Fig. [511 The extraploation results at k = are 
still close to the gs SCUT result and the extrapolation 
results at k = 0.5tt are in good agreement with the gs 
SCUT result. The linear extrapolation is even in ex- 
cellent agreement with our result. Note that the points 



from the L = 6 calculation are omitted for the extrap- 
olations because they deviate from the behavior of the 
other points due to the very small size. 

For x = 1 and A = 2 the SCUT result for the dis- 
persions comes close to the extrapolation of exact diag- 
onalization results in the regions where the continuum 
is distinguishable from the hole dispersion. This is most 
obvious for the odd band at k = 0. The left panel of 
Fig. [211 shows the finite size scaling for this case using 
the same extrapolations like for x = 1 and A = 1. We 
determine £ to be w 4.03 by the extrapolation based on 
exponential saturation. The correlation length is again 
used to apply this kind of extrapolation also for k = 0.5tt. 

The finite size scaling for k = 0.5ir in case of the odd 
band is investigated by both extrapolations in the right 
panel of Fig. [55] For both momenta the linear extrapola- 
tion comes close to the SCUT result, but also the extrap- 
olation based on exponential saturation does not exhibit 
large deviations from the SCUT result. The deviation for 
k = is rj 2% and the deviation for k — 0.5tt is s» 7%. 
The points from the L = 6 calculation arc omitted for the 
extrapolations because they deviate from the behavior of 
the other points due to the small size of the system. 

For x = 1 and A = 3 the tendency of the finite size 
scaling at k = 0.5tt indicates that the exact diagonaliza- 
tion results overestimate the energy of the odd band in 
this region so that one can expect that a proper finite 
size scaling yields a dispersion that lies completely below 
the upper boundary provided by the restricted SCUT 
energies. 

However, an extrapolation from the two points (L = 8 
and L = 12) at k = 0.5tt is not conclusive. We have 
several points for an extrapolation for k = and k = w, 
but there the distinction between continuum and disper- 
sion is difficult. Because the gs SCUT diverges without 
restriction of the generator, we actually expect that a 
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strong overlap of the energies of single holes and of the energies of hole plus triplon states is present. 



1 E. Dagotto and A. Moreo, Phys. Rev. B 38, 5087 (1988). 

2 E. Dagotto, J. Riera, and D. Scalapino, Phys. Rev. B 45, 
5744 (1992). 

3 T. Barnes, E. Dagotto, J. Riera, and E. S. Swanson, Phys. 
Rev. B 47, 3196 (1993). 

4 E. Dagotto and T. M. Rice, Science 271, 618 (1996). 

5 M. Troyer, H. Tsunetsugu, and T. M. Rice, Phys. Rev. B 
53, 251 (1996). 

6 S. R. White and D. J. Scalapino, Phys. Rev. B 55, 6504 
(1997). 

7 E. Jeckelmann, D. J. Scalapino, and S. R. White, Phys. 
Rev. B 58, 9492 (1998). 

8 J. Oitmaa, C. J. Hamer, and Z. Weihong, Phys. Rev. B 
60, 16364 (1999). 

9 O. P. Sushkov, Phys. Rev. B 60, 14517 (1999). 

10 D. C. Johnston, M. Troyer, S. Miyahara, D. Lidsky, 
K. Ueda, M. Azuma, Z. Hiroi, M. Takano, M. Isobe, 
Y. Ueda, et al, condmat/0001147 (2000). 

11 M. Brunner, S. Capponi, F. F. Assaad, and A. Muramatsu, 
Phys. Rev. B 63, 180511 (2001). 

12 C. Jurecka and W. Brenig, Phys. Rev. B 63, 94409 (2001). 

13 C. Jurecka and W. Brenig, J. Low Temp. Phys. 126, 1165 

(2002) . 

14 K. P. Schmidt and G. S. Uhrig, Mod. Phys. Lett. B 19, 
1179 (2005). 

15 G. Roux, S. R. White, S. Capponi, A. Lauchli, and D. Poil- 
blanc, Phys. Rev. B 72, 014523 (2005). 

16 P. Bouillot, C. Kollath, A. M. Lauchli, M. Zvonarev, 
B. Thielemann, C. Riiegg, E. Orignac, R. Citro, 
M. Klanjsek, C. Berthier, et al., Phys. Rev. B 83, 054407 
(2011). 

17 A. B. Harris and R. V. Lange, Phys. Rev. 157, 295 (1967). 

18 F. C. Zhang and T. M. Rice, Phys. Rev. B 37, 3759 (1988). 

19 K. P. Schmidt and G. S. Uhrig, Phys. Rev. Lett. 90, 227204 

(2003) . 

20 C. Knetter, K. P. Schmidt, and G. S. Uhrig, J. Phys. A: 
Math. Gen. 36, 7889 (2003). 

21 C. H. J. Oitmaa and W. Zheng, Series Expansion Methods 
for Strongly Interacting Lattice Models (Cambridge Uni- 
versity Press, Cambridge, UK, 2006). 

22 T. Fischer, S. Duffe, and G. S. Uhrig, New J. Phys. 10, 
033048 (2010). 

23 T. Vuletic, B. Korin-Hamzic, T. Ivek, S. Tomic, B. Gor- 
shunov, M. Dressel, and J. Akimitsu, Phys. Rep. 428, 169 

(2006) . 

24 M. Uehara, T. Nagata, J. Akimitsu, H. Takahashi, N. Mori, 
and K. Kinoshita, J. Phys. Soc. Jpn. 65, 2764 (1996). 

25 S. Gopalan, T. M. Rice, and M. Sigrist, Phys. Rev. B 49, 
8901 (1994). 

26 M. Matsuda, K. Katsumata, R. S. Eccleston, S. Brehmer, 
and H.-J. Mikeska, Phys. Rev. B 62, 8903 (2000). 

27 K. P. Schmidt and G. S. Uhrig, Phys. Rev. B 75, 224414 

(2007) . 

28 P. W. Anderson, Phys. Rev. 79, 350 (1950). 



29 S. A. Hamerla, S. Duffe, and G. S. Uhrig, Phys. Rev. B 
82, 235117 (2010). 

30 R. Coldea, S. M. Hayden, G. Aeppli, T. G. Perring, C. D. 
Frost, T. E. Mason, S. W. Cheong, and Z. Fisk, Phys. Rev. 
Lett. 86, 5377 (2001). 

31 Y. Mizuno, T. Tohyama, and S. Maekawa, J. Low Temp. 
Phys. 117, 389 (1999). 

32 S. Brehmer, H. Mikeska, M. Miiller, N. Nagaosa, and 
S. Uchida, Phys. Rev. B 60, 329 (1999). 

33 S. Sachdev and R. N. Bhatt, Phys. Rev. B 41, 9323 (1990). 

34 S. Duffe, Effective Hamiltonians for Undoped and 
Hole-Doped Antiferromagnetic Spin-1/2 Ladders 
by Self-Similar Continuous Unitary Transforma- 
tions in Real Space (PhD Thesis, available at 
http : //tl .physik . tu-dortmund.de/uhrig/phd.html 
TU Dortmund, 2010). 

35 S. D. Glazek and K. G. Wilson, Phys. Rev. D 48, 5863 
(1993). 

36 F. J. Wegner, Ann. Physik 3, 77 (1994). 

37 S. Kehrein, The Flow Equatin Approach to Many-Particle 
Systems, vol. 217 of Springer Tracts in Modern Physics 
(Springer, Berlin, 2006). 

38 A. Mielke, Ann. Physik 6, 215 (1997). 

39 A. Reischl, E. Miiller-Hartmann, and G. S. Uhrig, Phys. 
Rev. B 70, 245124 (2004). 

40 G. S. Uhrig and H. J. Schulz, Phys. Rev. B 54, 9624(R) 
(1996). 

41 S. Trebst, H. Monien, C. J. Hamer, Z. Weihong, and 
R. Singh, Phys. Rev. Lett. 85, 4373 (2000). 

42 C. Knetter, K. P. Schmidt, M. Griininger, and G. S. Uhrig, 
Phys. Rev. Lett. 87, 167204 (2001). 

43 A. Mielke, Eur. Phys. J. B 5, 605 (1998). 

44 G. S. Uhrig and B. Normand, Phys. Rev. B 58, R14705 
(1998). 

45 C. Knetter and G. S. Uhrig, Eur. Phys. J. B 13, 209 (2000). 

46 S. Dusuel and G. S. Uhrig, J. Phys. A: Math. Gen. 37, 
9275 (2004). 

47 T. S. Nunner, P. Brune, T. Kopp, M. Windt, and 
M. Griininger, Phys. Rev. B 66, 180404(R) (2002). 

48 The conservation of the parity with respect to the sym- 
metry axis of the ladder forbids the change of the triplon 
number by an odd value as the parity of one triplon is odd. 

49 The boundaries of this continuum are constant because the 
triplon dispersion is also constant for x = xu = 0. 

50 Even if the ROD exibits a small hump at small values of I 
before it decreases with a constant rate, the sorting usually 
does not pose a problem. 

51 Considering the momenta where the triplon dispersion and 
the triplon-hole-continuum can be clearly distinguished, 
the deviations between SCUT and exact diagonalisation 
are diminishing even further with growing system size, see 
the finite size scaling in Appendix IA1 



